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Nous prposons une construction hypothetico - deductive de la transformation 
de Lorentz dans le cas general et une explicitation des 2 termes de la forme 


polaire de cette transformation . 


Chapitre I : Axiomatique + Matrices de Lorentz ° 


On considere notre espace physique assimilé à un espace affine réel à 3 dimensions E 

dans lequel peuvent se mouvoir librement par rapport à E et entre eux 

des points appelés "observateurs " .On considère P un ensemble de ces points O 
 — — 


auxquels on associe un repère orthorme Ro| O, io jo: ko) etune horloge H,,. 
On suppose que les hypothèses suivantes sont vérifiées sur P : 


(Hypothese 1) : Pour tout point O € P, on suppose que pour tout point Q défini 
à partir de R,, on peut lui associer une horloge synchronisée avec celle de © . 


Conséquence : La synchronisation des horloges dans chaque repère permet d'associer 
a chaque observateur O un repère R, quadrimensionel d'origine © : O au temps t,=0 


puisque la variable temporelle t est indépendante des variables spatiales . 
On définit ainsi une famille de repères quadridimensionnels (Re) < 


R, = Ro(0 pre =C°t,e,=ke,;=J, a+) avec c la vitesse de la lumiere, 
0 


T vecteur temporel unitaire porté par la droite d'univers associées à ©, 
 — — 

et1, J, K sonti,, jo; ko au temps t,=0. 

On définit aussi sur R, la forme bilineaire symétrique définie par 

p(e; €) = pouri Æj, P(e e;) = 1, pe, €) =- 1 pour i=2, 3, 4. 

(Hypothese 2) : Si une ligne d'univers est une droite affine pour R, , 


elle est aussi une droite affine pour 
tout autre Rey . 


(Hypothese 3) : 
On suppose que les espaces E, défini par(R, ‘ H5) sont munis des mêmes lois physiques 
qui permettent de définir les mêmes unités de mesure d'espace et de temps pour chaque E, . 


On choisitun point arbitraire O,, de P et le temps oo, = ( indique par Ho : 


Pour tout point O de D, par translation spatiale prenons comme nouvelle origine 
spatiale de E, le point O'de R,, qui coincide avec O, au tempst, =0. 
0 


Prenons comme nouvelle origine temporelle de E,, l'instant ou O'et O, coïncident. 

Par la suite, on supposera que tout les repères spatiaux R,, auront leur origines confondues 
au temps t = 0 de leur horloge respective. 

Et par conséquent la famille de repères quadridimensionnels (Re) aura la même origine 


spatio — temporelle. 


(Hypothese 4) : On suppose que la lumiere a un mouvement rectiligne uniforme dans tout E,, 


et que sa vitesse numérique c est constante et indépendante de la source d'émission 
et de l'espace E,, considere. 


Conséquence : Si X = (ct, x, y, z) est le vecteur représentant la position 
d'un photon issu de l'observateur O au temps t = 0 pour le repère R, 


etsi X'= (ct',x',y',z') est le vecteur représentant 


la position de ce même photon issu de l'observateur O'au temps t'= 0 
RE TN 0 
ona = À > LEA RE pourt >0ett'>0. 
t' 


Donc les coordonnées du photon vérifieront simultanément : 
PRET SR ES DD De 2 
CC-X HY + =0 et CE —x" +y" +z “0dans RetR,, 





Si on considere la forme quadratique associée à @ : 
DIX) = CP -X +ÿ +: DX) =0 = D(X') =0. 


(Hypothese 5) : Isotropie de l'espace . 
Nous énoncerons cette hypothèse au chapitre II. 
( Voir auussi C. Semay , B. Silvestre — Brac : relativité restreinte p.105) 


Le probleme à resoudre : 
On cherche une transformation T : R4 — R4 inversible qui donne les coordonnées spatio 
— temporelles 
du point P pour l'observateur © en fonction de celles de l'observateur O", 
telle que sous les hypothèses I, 2, 3, 4 ci-dessus possède les invariants suivant : 
(1)T(0)=0 , 
(2) T transforme les droites affines en droites affines , 
(3)si X =T(X") : BD(X) =0 = B(X') =0 c'est a dire que T laisse invariant le cône d'isotropie de ®. 


Remarques : 

L'inversibilite permet d'avoir un groupe de transformation . 

La premiere condition vient de l'origine commune des 2 repères et que des que T est affine, 
elle est linéaire . 

Les 2 autres conditions traduisent les 2 invariances admises par les hypothèses 2 et 4. 


4 
Le seul paramètre qui lie O et O'est la vitesse relative V de O' par rapport a O et donc T = T( V). 


Pour cela on va utiliser les 2 théorèmes suivants donnés dans toute leur généralite : 

Théorème 1 : ( P + Boyer ." Algèbre et Géométries ". C&M 2015 p.51). 

Soit f: 6 — Eune bijection d'un espace affine & de dimension > 2 sur un corps K 
possédant au moins 3 éléments, telle que pour tout a, b, c alignes alors f(a) , f(b) f(c) sont alignés, 
alors f'est semi — affine. 

SiK=R fest afjine. 

(Ce théorème est appelé theorème fondamental de la géométrie affine + Abouabdillah Driss a proposé 
une version améliorée de ce théorème : (https //ssrn -com/abstract = 3181422) en remplaçant 
l'hypothèse f bijective par f surjective ). 

Conséquence : Dans un espace affine réel les bijections qui conservent les droites affines 

sont les applications affines bijectives . 

On note C($) = {x/@{x, x) =0} le cône d'isotropie de ® une forme bilinéaire symétrique 

etRad = {x/V y Q(x, y) =0). 

Remarque : Si @est la forme bilinéaire symétrique associée à D= cf -x +ÿ +2: Rad $= {0}. 

Théorème 2 : (R.Goblot ."Algébre linéaire "Masson 1995 p.254) 

Soit: E x E —K, E espace vectoriel sur K, une forme bilinéaire symétrique telle que Rad @ Z C() 


Pour qu'une forme bilinéaire @" soit proportionnelle à @ il faut et et il sufjit que C($) = C(Ÿ' ). 


(Comparer ce théorème associe aux Jormes linéaires ayant un noyau commun .) 
Les conditions (1) et (2)imposees à aT implique que T est linéaire ( Théorème 1) . 
Soit M est la matrice qui représente T dans les bases associées aux repères R, et Re 


Sionnote X'=MX ,ona XX) =-0 = BX"') = BMX) =0 
alors (Théorème2) :11 70, VX D(MX) =1- B(X) 


110 0 0 
0 -1 0 0 
en notant G= On a : 
0 0 -1 0 
0 0 0 -1 


{MX)G(MX) = 'X MGMX =2-XGX = MGM = 2-G 


la dernière équivalence est conséquence de la symétrie des matrices MGM et G. 
(J-M. Monier."Algèbre 1 et 2". Dunod 1997.) 


Evaluation de À : 


— 
Montrons d'abord que si la vitesse de O' mesurée par observateur OestV, 


la vitesse de O mesurée par l'observateur O'est -V: 
(Pour le moment on ne dispose que des hypothèses ci — dessus ) 
(cf C * Semay, B -Silvestre - Brac + "Relativité restreinte". Dunod 2010 p108 notes ) 
On se place d'abord du point de vue de l'observateur ©. 
Considérons le cas ou O"s'éloignant de ©. 
Plaçons nous du point de vue de l'observateur ©. 
Un rayon lumineux part de O vers O'au temps t, , atteint O'au tempst, au point P,, 


et revient immédiatement vers O, qui est atteint au temps t, puis repart immédiatement 
vers O" qui est atteint au temps t,au point P,, 
puis retourne immédiatement vers O qui est atteint au temps y. 


Il y a donc un double aller — retour. 

La vitesse numérique c du rayon lumineux étant la même dans les 2 sens on a : 
2 ; 2 

Lorsque le rayon lumineux atteint, au tempst,, le point P,, 








by — 


O" continue a s'éloigner de O alors que le rayon retourne vers ©. 





A PU est la durée entre les 2 contacts du rayon lumineux et O'. 


> — ty = n 
Si on suppose que O" a une vitesse uniforme V de module V on a | PP, | =V 


— 
On peut remarquer que | PP, | est aussi égal à la la différence 
UD | 





— — 
or] - [or] -[ 2 - 

















1, —21t, +t 
Donc V =c : 2 “el - 2 L | 
ty Lit 
2 
” (4-6) —2(6-4) se mit 
l, “lo PE) ; 


Soit M = (mn; j) la matrice de transformation telle que X'= MX, 
X les cooordonnées d'un point P dans la base de R, : 


X"' les cooordonnees d'un même point P dans la base de Rey 


Plaçons nous maintenant du point de vue de l'observateur O" qui veut évaluer V 
Considérons les coordonnées de l'observateur O dans R, et Rey ; 


O a pour coordonnee (ct, 0, 0, 0) dans la base assciee à Ro ; 
et Oaura pour coordonnées : 

(a',x!,y!,z') = (mn, 1m cts My °C, My] *ct) pour la base assciée à R y 
Donc ct'=m 1,1 °Ct: on peut donc écrire qu'en O' t=@-t" avec nécessairement & > 0 . 


Si l'observateur O' évalue V'" en observant l'horloge sur l'observateur ©, 
il obtiendra la même valeur que l'observateur © puisque en faisant le même calcul il obtient : 


Qt, — Qt t, -t 
V'=cl 1 —2 2 espere, 
œ-t, — ty tt 








La vitesse dans les 2 cas est collinéaire à O0" et le sens de V est celle de O0" pour O 
et O'O pour 0" 


On peut maintenant évaluer À : 

(NM.J. Woodhouse." Special Relativity " - Springer 2002 p.80) 

D'apres l'axiome 2 les lois de la physique étant les mêmes pour les espaces associées 

à O et O" la dilatation des durée sera la même que ce soit l'observateur O observant O' 
ou bien que ce soit l'observateur O' observant O puisque les situations physiques sont les mêmes, 
le modules des vitesses étant identique. 


Les coordonnées spatio - temporelle de O0" sont représentée par le vecteur X = (ce, x, }, 2) 


pour l'observateur ©, et par X (ct, 0, 0, 0) pour l'observateur O". 
Comme X = MX" on en déduit que le temps t mesuré par un observateur situe en O d'une horloge 
située en O'et qui indique le temps t' à un observateur situe en O' verifie t=m, jt! 


avec m,, > 0 sion suppose qu'iln'y a pas de retournement du temps car : 


ct , 
_ L'EAU Dr) my, 1Ct 
B,ct m m m m m, ,ct! 
+ 2,1 "22 "23 "2,4 0 2,1 
y VC My Mj2 M;3z M4 0 m3 1 Cl 
z 


0 , 
Bet Mes MES ee MAeE 


Considérons la situation symétrique où un observateur situé en O' observe une horloge 
situee en O qui indique un temps t pour l'observateur situe en O -* L'observateur situé en O' 


mesure alors un temps t' - Comme X"' =M !X et sion poseN=M ! on aurat'=n; jt 
etn; , > 0 comme précédemment 

Comme nous sommes dans la même situation physique que précédemment n 11-11. 
De'MGM=À-G ,À Z0,onen déduit que M=À-G-M /G! 

D'où M = "(G'M-G) =m, ,=n,,=(M),,=A (G.M-6), ;. 

Or (G-"M-6), 1= M1 donc À=1. 


Les matrices M qui vérifient MGM = G sont appelées matrices de Lorentz 
et forment un sous — groupe L du groupe GL, (R) * En particulier M !=G-M:G et det”(M) =1. 


Pour expliciter M dans le cas qui nous intéresse on va utiliser les 2 théorèmes suivants : 
Théorème 3 : (J-M. Souriau."Calcul Linéaire " - PUF 1964 p.378, voir Annexe — ) 
Toute matrice M de Lorentz peut se mettre sous la forme : 


£ 0 
M=exp( aN) - oùaEeR,E£=+1, 
0 Q 
0 ; 3 t t 
N= ,avecX E R'telque : XX=1, Q00=Id,; . 
X 
ch( æ) sh(æœ)'X 
De pl aN) = 
CR (ax (4, +(eh( a) -1)XX) 


Théorème 4 : 
(R. Mneime , F. Testard "Introduction à la théorie des groupes de Lie classiques ". 
Hermann Paris 1986 + Voir annexe — ) 


Soit M une matrice inversible a coefficients réels alors il existe un couple unique de matrice 
S'' symétrique définie positive et O orthogonale telles que M=OS". 


Conséquences : 


t 
On en déduit en posant S=0S'"O < S'='OSO que M peut se décomposer de manière unique en 
M=0S"=0('0$0) =$0. 
Cela entraine que décomposition du théorème 3 est unique : 
il suffit de vérifier que exp( aN) est définie positive . 
oi est symétrique réelle , il existe donc une matrice U orthogonale réelle et une matrice diagonale 
réelle D = (d; ;) telle que aN = 'UDU. 
Comme exp(aN) = Uexp(D) U et exp(D) = (exp(d, i) ) les valeurs propres de exp( aN) sont 
strictement positive et exp( aN ) est définie positive . 


De la formule det( €) =e" 4) on en déduit que det(exp( aN) ) =1 
et que det(M) =æ-det(Q) =-+1. 


Chapitre IT : Expression de la partie symétrique. 


On considère maintenant un point O'qui s'éloigne d'un point O à la vitesse constante V.. 
— 


Onsoteh = " où c est la vitesse de la lumiere. 

On considère les 2 repères Re! R définis plus haut . 

X étant les coordonnées d'un point P dans la base associée à R, 

X ‘étant les coordonnées de ce même point P dans la base associée à Re : 
M la matrice définie par X = M ° X' la matrice de passage de R, d Re ; 
M est donc une matrice de Lorentz et s'écrit donc sous la forme 


€ 0 
M=exp( aN) - ouaŒER,E=+1,N-= 
0 €Q 


x 
X oO 








avecX € R° tel que : 'XX=I, ‘20-14, . (voir théorème 3). 


—_ 
On va exprimer dans un premier temps exprimer exp( aN ) en fonction de B. 


Théorème 5 : 
En reprenant les notations du théorème 3 on a, si y= : 
2 
1-B 
tr > 
Y 18 | 
(aN) > tt > avec LB | les coordonnées de P dans R 
exp = D 
hé) 0, + PÉThEl 
(1+7) 


Démonstration : On peut donc ecrire que les coordonnées de O'dans la base associée à 8 ,, sont : 
— 


t t t Re 
W='{ct,1V,1V,tV,) =ct (1 B, = — 

(ct, 1Vp1V}1V3) =c ( B, B, B;) avec B À 
et dans 8", :'W'= “(ct', 0, 0, 0) = ct'"(1, 0, 0, 0) avec W=M-W. 








1 0 
On remarque d'abord que W'=W'donc W=exp(aN)W". 
On est ramene a : 
1 
1 
ch( à) sh(æ)'X ||, B, 
— t , 
sh (æ)X (4, + (ch( æ) -1 XX) 0 B, 
0 





1 ch( à) 
B, sh (æ)X, 
cela entraine que t =1 et donct=ch( æ)t' d'où : 
2 sh (æ)X, 
B, sh (@)X, 
1 ch( @) 
Flo (a)B=sh(a) 
ch\ = pourt'£0 = ch{ @a)B=sh(a)X, 
B, sh (æ)X, B 
B, | |sh(a)x, 
> > 2 2 2 
donc B=th(a)X = B =th"(æ) puisque X =1. 
Remarque : 
— 
: 22 LÆ Es 
Comme |th(æ)| <1 V@œÆEëR, nécessairement Ip = | | <1< | v | < € 
ce qui exclut toute vitesse supraluminique entre observateurs . 
—2 1 
OnposeB=\ B et y- —— 
2 
1-B 
1 2 1 1 2 
Comme 1 — th°( @) = ——— = ch @) = = =Y, 
ch°( @) 1—th(@) 1-B 


comme ch( œ) >1 y=ch( @);commesh” (æ@) =ch°( @) _1=Y — 1 


nd mis F ph. 








—2 
1-B 1-F 
En résumé on a y= ch( @), YF = sh? (æ), p =1h(@) : 
D'autre part : 


2 
YB=Y —-1=(y+1)(y-1) = ARE = ch( @) —1 








(1 +7) 
__ (BB) _ (BB) 
et XX; = (a) = F donc 
2 BB. . 
(ch( @) —1)XX;= r p ) = y BB.. ce qui permet d'écrire que : 


(+7) (1+7y) 


A, +(ch( @) -1)XX 


2 
a er (ee rt 
G+pfr (Ginrr BB 






































(1+7) 
2 2 2 
losn# Gent Gent 
2 2 2 
Gen Gene Gan ls 
> tr > 
y. + LBlDB! 
3 
(1+7) 
sh (æ)B 
Comme sh (@)X,= = =ch( @)B.=7B. d'où le résultat . 
th(@) ne 
En résume : 
tr > tr > 
T 118 | L 0 y [16 le 
M = > tr > = > jte > 
| nu, + BI 00 LB] nu, + Bo 
(1+7) (1+7) 
avec ‘20-14, . 
h8 l'I$] y él 
Par la suite si on note C=Id,z + — M = F 
(1 +7) LB | 26 
Remarques : 
_ > {6 > Le & > >> " 2 
oclfl=lr, + Lélhél ÿ-6 +7 l8LIAIBI [Gifs 4 Ÿ | 
U+r) (1+7) (1 +7) 
fre ÉL) rl 
(1+7) 
(2)Si M=M"M" est le produit de 2 matrices de Lorentz sans rotation avec : 
t, — t, — t, —> 
nn (6 )e el, & (yB") DAT (y"B") 
> — — 
Y ca yB' C:' vy'B" GC" 




















> — > t—> 
y v"CB'B" +1) y Y'B" +yB'C" 8 > 
alors M= = y-yy"(B'B" +1), 
— — 6 
YY'B' + y'C' B' y Y'B' B r + C ’ C" 


B'+c'p" 66 

B' +C'B" : S 

TE .C=1d + - a a=C (y yB'B"+0c'c"). 
y (86"+1) 





- 


le bre dre 





61116 BB 
Id. + 1 El 
F (1+7) Æ  (1+7) 














(147) (1+7) (147) 
BB 
2 2 
= gr +7) +7 +7) 77 —1)) =. 
(1 +7) 
-1 
tr > 8 
Donc C'} = 3 + él b8l =1d,, - ve» cequipermet d'écrire Q 
> > 
en fonction de B', B'" : 
6e 
o=|1,- y PB (y yBrB" +c'c") avec: 
(1 +7) 
— B" +c'B" 2B'B np" p" PB 
B= ,C'=1 ++ ac"=id + TRE y y y BB" +1) 
8 
(3)La connaissance de C'° =Id,, - Y.——— permet d'évaluer Qen fonction de M : 
de C'=ld,z FRS d'évaluer Qen f d 
y 


Méjoy My is Mig Mrs 
MES Mrs" MS Z 
Mg Miss Ms M3 à 
M= + Considerons le bloc M=| M; M;z M, 
MSc. Miss Ms:e SES 
MS y Mig Ms 
My Ms Mis Mis 


qui représente les composantes spatiales des vecteurs spatiaux de la base associée à O' exprimées 


dans la base associée à O. 


Onaura A=C M. 


> 
4)Si B/ i alors les termes non diagonaux de C sont nuls , 








; 2 (14+Y)+%8) 
etun seul terme diagonal de C est different de I : 1 + : F = 








(1 +7) (1 +7) 
= ER ER ne = car ÿ 1 = L 1 = b -ÿ$. 
(1 + ) 2 —2 
: 1 —B 1—B 


2 2 1 
(5)Comme y F =y —1on peutremplacer , en notant 8, le symbole de Kronecker, 


] Ÿ BB; rs ; t 
DE G " = BB; pa ar & + (y — 1) F dans l'évaluation de Id: + (ch( @) -1)X'X 











e Fa 14 _B : El 7 -1 . 


et on retrouve que argth( B) =argcosh(y) =æ car y=ch( @) etB=th(a@). 


(7) Sachant que |exp( aN) ] "=exp( (-a)N),7y=ch( a@)et B=t(a)X il vient : 


y hf] U | 
: > tr > = ss > tr > 
PL He + (an LS TT 


(8) Valeurs propres d'une matrice de Lorentz symétrique : 
Tout calcul fait on a : 


















































RES M D... 
7 B 00 2 2 7-1 +B) 0 0 0 2 2 
0 0 0 1- 0 0 
AP 7 dE A se) PE NE, SE, y 
0 9 0 1 0 2 2 0 0 1 0 2 2 
0 0 0 1 0 0 10 0 0 0 1 0 0 1 0 
0 0 0 1 0 0 0 1 
X Æ LA LA 
+) sn( +) 0 0 y +B) 0 0 0 +) sin( +) 0 0 
+ Es 0 7U-B) 0 0 x x 
=| sinl cos 7 0 0 | , à vs: sn) +) 0 0 |; 
0 1 0 0 hi 0 0 0 1 0 
0 0 1 0 0 0 1 
Donc les valeurs propres de À, sont y-(1 + B), y-(1-B),1. 
En se servant des remarques du théorème 5 on a : 
+B 1-B 
-LTHD)= et donc que y-(1-B) = à 
7 (1 +B) TB que y-(1-B) 1+B 
1 1 +B — 
De plus EL 1=B = argth(B) =@, les valeurs propres de M" peuvent s'écrire : 
1 1 1- 1 : 
+B EXD ne Lb Sea, b =eXp RE ER =e et 1. 
1 —B 2 (1-8 1+B 2 |1-B 
Plus généralement : 
tr > 
Y 18 | 
Calcul des valeurs et vecteurs propres de M = co > tr > 
Hô] nu, + 6Tbél 
(1 +7) 
On remarque que : 
tr > 
u CE B+Y 
Dir > > = LB tr > 1 
à 21.181. LB1l# 18] 
FI (BIT | #6]+[8]+ 
#  (1+7) (1+7) 






































se FEAR. LS | 
pl +rlél | | U+B 8] "AI 
7 #6] _p -B+% 
de même : > ir > = tr > 
B 21.131, 16116116] 
né] ns + DÉLOÉL LIBIT Lg] + [814 HÉCDÈDE 
B+ | B+h 
CPE OS EETIPTE PET 
OR PRES j 
Bhél+ 716] | | U-BrËl HDI 

















Le produit des 2 premières valeurs propres est P'= (1 +B)y(1-B)y=1, leur somme 
est S'= 27 * Or det(M) = 1 et Trace(M) =y+3 +y-1=2(y+1). 

Donc la somme des 2 autres est S = 2. Donc les autres racines sont égale à I. 

Si on a : 




















tr > 
Y 8 | 
a a 
> tr > à] = à] alors 
6] 4 ' d 
F" (AP) 
tt > 
m+ [#14] 
a 
D it > —|| 
RL [a] [a] 
al] +141+ 
(1 +7) 

Une solution possible est a =0 et tout vecteur u tel que [BITu] =. 
0 0 

2 vecteurs propres possibles associés a la valeur propre 1 sont |  |et | , 
u y 

tr tr 0 
avec [Bllu]=0, [Bllv1=0, vl[u]=0. = let | = | sont du genre espace. 
u v 














—+ — 
Plus generalement tout vecteur de la forme | _, | avec w L B est invariant. 
w 


Lé 
On aura l'ensemble spectral suivant associe a M = S tr > 
— 
Hé] 1, + 161 
(1+7) 

B . 0 0 

> , (1 +B)7 |, > , (1-B)y|, > LT}, > L1 : 

[8] [8] u v 














Remarque : 
On démontre plus généralement que toute matrice de Lorentz M a au — moins 


une valeur propre @ réelle associée a une direction propre isotrope 
(E. Gourgoulhon , relativité restreinte, EDP Sciences p.183). 
Le polynôme caractéristique x(À) de M se factorise en (À — æ@)Q(À) où Q 


est un polynome du 3°"* degré et donc M a au moins 2 valeurs propres réelles . 


1 0 0 0 
0 1 0 0 
La matrice de Lorentz : . montre que 2 est une valeur maximale dans le cas 
0 0 cos(@) -sin(@) 
0 0 sin(@) cos(®) 


général. 
Toujours dans le cas général, si une valeur propre À est associée à un vecteur propre non isotrope X, 
alors À =+ 1 car si MX =AX = 'X'M=XX on a, comme ‘X 'MGMX ='X GX car M est de Lorentz, 


et on a aussi IX MGMX =X (X GX). 


(9) Autres formulations : 
On a vu qu'une matrice de Lorentz M sans rotation pouvait s'ecrire sous la forme 
_. 


| 








Lé 
M = > tr > c'est a dire sous forme développée : 
6] 2, + 611 
(1+7) 
tr > 

7 Y [6 | ax 
X=MX'S ee = > sfr > > 

L Hé] 2, + 6161 || [rl 











On trouve dans la littérature d'autres formes par exemple : 




















(a) forme symétrique : 
# Le tr — 
[cr] = | cr] +[# | ct'+|r ] 
(1 +7) 
> > re > 
[r]={r"] +11] [Le] + cr’ 
(1 +7) 
Pour cela on remarque que : 
tr >> tr > 
a=yt'+ [8 Îlr'l=cr"+ no 198 ÎLr'l or 
2 2 y 
comme ÿ ff =Y —1  ona y 1: 
(1 +7) 


(b) forme avec double produit vectoriel note x : 


a=yer+ TBllr) 


= y([Bler+ L)) + 2 (Rx ([B1x[r1)) 








On remarque que : (= = 


[B| x ( [8] x pal = )B | BI). -( formule du double produit vectoriel) 


et que 


AS Ce El ANERE eue a 
a 


y (TB ler+ [r°]) + (Rx [8x 1) 
de Me 2 le). 81 (TRI). 1) 


(y) - CD TG$] 








= yBler + [1 + [Fr] 











F y 
ile ie BB) 246 RÉ 4ee)- cm 








Chapitre III : Expression de la partie orthogonale . 


Maintenant éevaluons £et Q 


Evaluation de £=+# 1 : 


Comme 
W=MW "= é [2 À ' W'= . LB le W" 
16 | C 6 el | ca 


avec W='(cr, (VV, 1V;) , W'= l(ct!, 0, 0,0) .Onen deduitque t=£-yt'. 
Donc si £=-1, il y aurait un renversement du temps difficile à justifier physiquement. 
Par la suite on suppose que £=+ 1. 


1 0 
0 Q 


A 
Reste à évaluer @ . On posera Q = 








Remarquons tout d'abord que M est symétrique = Q = Id}; 


A 
carM=SQ =M Id, avecS symétrique définie positive , par unicité de la décomposition Q = Id}. 


Evaluons d'abord directement un cas particulier : 


Theoreme 6 : 
SAS OT 


Ye, => 
Dans lecasou V/Ox// 0% x'et de mêémesens, 0/0! O'y'et de memesens, O07z/0"z' et de mêmesens, 


7 YB00 
_ er d 
M se réduit à À,= FT 04 avec B = >» B=V B-B ,7- 1 


0 0 10 € Fer 


0 0 0 1 
Démonstration : 
Le _ 

Si DAOx 0%" de mêmesens, 0/0! ‘de mêmesens,, O0z/0'z' de mêmesens, 

7 Y7B 00 
A : 

M s'écrit sous la forme M=A,.Q avec À ,= TB 7 00 . 

0 0 1 0 


0 0 01 




















A 1 0 
et Q= où Q est une matrice orthogonale d'ordre 3, 
” 7?2 2 
PS je (z +y)+78) L'ÉVEV EL. 
(1+7) * (1 +7) (1 +7) 
B__r? 
2 1 
cary —1= nd = =7Y/B. 
1 —B 1 —B 
1 0 0 0 
0 0 
A 0 ®,, D; ®,, 1 0 
SiX=M -X"'etsi Q = Xi » À'}= 
0 ©@,, @,, @,, 0 1 
0 0 
0 ®,, @,, ®,, 
YB@, , YB©, ; 7B@, , 
A Fe A re A V® 4 
A,9X = ; A,9X 5 et A,9X = 
®, ) ®, ; ®, y 
@®, 2 os ®,4 
Y®, ; 1 FR: 0 
Les hypothèses impliquent que | @,, [7/1 0 |,| @, |/\1|e 
0 0 
es ne 
1 0 0 0 
k ®, ) 0 0 
d'ou : Q = 
0 0 ®, ; 0 
0 0 0 ®, , 
0 Ra nr 
Comme @ est orthogona D, 70,3, @,,1 


et X", = 


nn OS S S 


on a : 


7 +yB 0 0 
A + 0 0 A 
D'ou A,Q SE B +7 les sens des axes étant conserves nécessairement © = Id, 
0 +1 0 
0 0 +1 


Bases standards et bases "parallèles" 

Considérons les points O et O'tel que O'estun mouvement rectiligne uniforme de vitesse V dans F. 
Munissons O d'une base orthonormale directe g(i i,j,k x) telle que i / Ÿ 00! , 

Considerons le plan P defini par ©, i et j. es Les Re 

Munissons O' d'une base orthonormale directe B (à ,J, k ) telle quei'/ V 00" etj' dans le plan P 
de même direction que j . 

Hypothese : Isotropie de l'espace. 

(Voir aussi C. Semay , B. Silvestre — Brac : relativité restreinte p.105) 

Nous traduisons l'hypothèse énoncée dans l'ouvrage cite ci — dessus , en l'éenonçant sous forme 


matricielle aux bases B et B'. 
Considérons une rotation d'axe OO" définie dans la base B : 


10 0 0 1000 

0 1 0 0 : 0 1 0 0 
R(@)= | etune symetrie S = 

0 O0 cos(@) -sin(@) 0 0 1 0 

0 O0 sin(@) cos(o) 0 0 0 -I1 


par rapport au plan Oxy (changement de l'orientation de l'axe Oz). 
Hip Es Mis Mr 4 


LT ST RE ne 
Appelons M = la matrice de Lorentz associee à O et O0": X = MX", 


mes 7, Mo: JS MES 2 


Var 4e Ps is 
L'hypothèse de isotropie de l'espace, dans notre cas particulier, se traduit par : 
R(e)-X=R(p) -MX'=M-R(p)X' ie R(p)M=M-R(p) oùM=R(e) -M-R(p) Ve. 
etpar SX=SMX'"=M - SX" ie SM=M°-S ouM=S°M°S. 
Si on fait une rotation d'axe OO" d'angle og =x: 


10 0 0 Far 42 25 ue | | 0 0 0 
01 0 0 Mot 22083 Mer 01 0 0 
0 0 cos(x) -sin(x) M; My Myz May 0 0 cos(x) sin(x) 
0 0 sin(x) cos(x) My My My My 4 0 0 -sin(x) cos(x) 


Mig Mg TMps My LS Rd ee EE 


Mpr. Ms. sé M5 M; gp Mi. M; My 
Mg ss Mise “Mer LUS ed LE be 
Ré ge ge Hs, MES My] My) M3 My, 
My 7 My 0 0 
M; 7 M) 0 0 
Donc M = 
0 0 M; 3 M4 
0 0 my 3 My 4 


On a en faisant maintenant une rotation d'axe OO" d'angle g =x: 


1 0 0 0 ; 
m m 
0 1 0 0 PT NE 
0 0 
0 0 cos à) sin[X ) sn. 
2 2 0 0 M; 3; M4 
0 0 sin( à | cos X) (à 0 Macs res 
My Mi 0 0 My My 0 0 
Mi, 7 M) 0 0 M, 3 M5 0 0 
0 0 My My 3 0 0 M 3 M4 
0 0 -m;, m;3 7 
d'où My 4 My 3 My y My 3 et donc : 
My M2 0 0 
M 1 M) 0 0 
M = ; 
0 0 M; z My 
0 0 M3 4 M3 
En apliquant l'égalité: M=S°M°S, 
0 0 
100 0 || "11 "2 100 0 
0 1 0 0 My M2 0 0 0 1 0 0 
0 01 0 0 0 m,, m,,||0 01 0 
0 0 0 0 0 0 -1 


É. 0 0 -m;, m3 


1 0 0 
0 1 0 
0 0 cos . 
0 0 sin 


M; 7 M) 0 0 
= et donc m; ,=0 et, 
0 0 m - mn 2 
3, 3 3, 4 
0 0 M; y M3 
M; 1 M2 0 0 
M; 1 M) 0 0 


> > > > > — : ; 
B({ i,j,k ) etB (ir, j", k ) ont donc leur composantes paralleles 2 a 2 pour les 2 observateurs O et O 
’ 


On peut donc appliquer le théorème 6 : 


7 YyB00 
M=A,- Y7B 7 00 
0 0 10 
0 0 01 


On appellera par la suite bases standards associées au mouvement relatif de O et O'Ies bases B et B' 
définies ci-dessus. 


On peut maintenant introduire la définition suivante : 


Définition : 2 bases B = (vip ki) et B'= (i 15 K; À associées à O et O", sont dites 
(à composantes) parallèles si le même opérateur orthogonal Q transforme simultanément B et B' 
en 2 bases standards associées à O et O'. 


Théorème 7 : 


Dans le cas où les 2 bases spatiales associées à O et 0": B et B'sont à composantes parallèles 2 à 2 
— 


> A _ 
(ici V n'est plus nécessairement parallèle à Ox) ,ona Q= Id}, et A(B) =RA, R avec 


1 


= à ou Q est l'opérateur qui transforme simultanément B et B'en 2 bases standards 


associées a OetO'. 

Démonstration : 

Notons 8 et B' les bases spatio — temporelles associées aux 2 bases spatiales à composantes parallèles, 
8, et B, les bases spatio — temporelles associées aux 2 bases standards . 


M la matrice de passage de la base 8 a la base 8". 


On a le schéma suivant: B —— @' à 
R | LR 


À, 


B — 7, 
1 
; ’ >, A 
CommeM=RA,R < À,= RMR et comme M est de Lorentz, elle est de la forme M = A( B) Q. 
t >, A > A 
En lisant le schéma ci — dessus À, ='RMR ='R (A B) O)R = (RA(B)R) ( ROR). 
> A 

On remarque que "RA B)r est symétrique et que 'R@R est orthogonal . 

t > A 
Comme À,= A ,: Id}, = ( RA( B)r) ( ‘RQOR ) , par unicite de la décomposition 


A  . A > ; 
ROR =Id,, . On en deduitque Q@ =1Id}, et que A(B) =RA,R. 


Comme 
| 1 0 0 0 7 7B00 1 0 0 0 
0 1 0 0 yB 7 00 0 1 0 0 
0 O0 cos(@) -sin(®) 0 0 10 0 O cos(@) -sin(®) 
0 0 sin(g) cos() 0 0 0 11||0 0 sin(g) cos(p) 
7 Y7B00 
TB 7 00 
0 0 10 
0 0 01 


A, est invariant par une translation autour de OO" et donc RA, R ne dépend pas 


du couple de bases choisies. 


Evaluation de Q dans le cas genéral : 


Theoreme 8 : 


Considérons maintenant O et O' munis de 2 bases spatiales quelconques B et B', 
Munissons O et O' de 2 bases standards B, et B",. 


Onnote B, B", 8 , B", les bases spatio — temporelles correspondantes 


SoitR etT les 2 opérateurs qui permettent de RACE respectivement de B a 8" 
— 


et de 8 à 7, * Alors A(B) = (RA, Te Q=R'T. 


Demonstration : 


On a le schema suivant : 


@ M, @ 


R J IT 


À, 


B —> à, 
1 

; ; >, A 

Comme M=RA,T < A,= RMT et comme M est de Lorentz, elle est de la forme M = A(B) Q. 
>, A > A 

En lisant le schéma ci — dessus A,= ‘RMT ='R(A(B) O)T= (RA(B)R) ( RAT). 
Comme À, = À, - Id _(RA(B)r) (RQT) cité de la dé ti 
omme À,=A,: Id,,= B , par unicite de la décomposition 
t A r 2 A t ne t 
ROT=Id, . Onen déduitque Q=RT et que A(B) =RA,R. 


Comme au théorème précédent R'T et RA, ÎR ne dépendent pas 


du couple de bases standards choisies. 


A > 
Conséquence: Q=M ( 0 ) 


Les opérateurs R etT sont indépendants du module de la vitesse relative de O et O'et 
Q est donc indépendants du module de la vitesse relative de O et O'. 


EN —> A ‘ 
Comme lim  A(B)=H, ona lim  M(B)-Q-=RT. 
IB| — 0 | IB — 0 
Si la vitesse est nulle O = O"indefiniment. 
On peut donc interpréter @ comme la matrice de passage de B à B' si la vitesse relative est nulle. 


Remarque : N °M + J + Woodhouse dans Special relativity p.82 propose la decomposition : 


n ga 1 0 1 0 
M=H.A,-K = A; lou HetK sont 2 matrices orthogonales 3 X3. 
0 H 0 K 








On peut alors écrire : M=HA,-K = aa x) (aix) =A-Q, 


par unicité À = HA, “À et HK=Q. 


Chapitre IV : 


#Exemples numériques : 


#Pour illustrer cette étude des matrices de Lorentz on va donner plusiurs exemples numériques 

#qui permettrons de vérifier la justesse des formules données précédemment . 

#Pour cela on va utiliser le logiciel Maple de Maplesoft approprié aux matrices de petites dimensions . 
#Liste des exemples : 


#(1)un proramme qui génére de façon aléatoires des matrices de Lorentz. 
#(2)Maintenant un programme qui génére 2 matrices de Lorentz symétriques et leur produit qui n'est pas symétrique. 
#(3)Décomposition polaire : M = A-Q.. 


#(4)Un programme qui donne les parties symétriques et orthogonale d'une matrice de Lorentz 
#par la formule :C-1.771 





t 
#(5)Transformation d'une matrice de Lorentz symétrique M sous la forme M = À + Ao.A 


#et d'une matrice de Lorentz M non symétrique sous la forme M = A :* Ao : B. 


1 0 0 1 0 0 
0 1 0 ns 0 1 0 e 
Sp. i . .B° BXT = a .B + BKT. 
#(6)vérification de la formule : + 1% B-B TE BB’ 
0 0 1 0 0 1 
1 0 0 
re . 0 1 0 RP 
#(7) Vérification de + 137 -BeBKT| .B=yR 
0 0 1 


#(8) Vérification de la formule (cas du produit de 2 matrices sans rotation ) : 


= RIT | aprem cr. 


#Exemples : 


#(1)un proramme qui génére de façon aléatoires des matrices de Lorentz : 
#Cas général : 


= 
#On génére aléatoirement le vecteur B ce qui permet de construire la partie symétrique de la matrice. 
SF — 


#On génére aléatoirement trois vecteurs ©, @2 et «3 etpar la méthode de Gram — Schmidt 


#on obtient un système ortonormal de 3 vecteurs ce qui permet d'avoir la partie ortnonormale 
#de la matrice. 


#En retour on aura la matrice gérale , vérification qu'elle est bien de Lorentz, la partie symétrique, 
#la partie orthogonale et la vérification que cette dernière est bien orthogonale et finalement B 
#On se limite à 10 chiffres significatifs. 


restart : 
with(RandomTools) : with(LinearAlgebra) : 


randomize( ) ; 


pl =pro( ) 

localG, co, B 3 M, A, Ex, By BR, ©, GS, ol, 2, 8 : 

unprotect(}) : 
1 0 0 0 
0 -I1 0 0 

G := 

0 0 -1 0 
0 0 0 -1 


co := evalf(\[ 7) : 


Generatel float(range =-1 ..1, digits = 10) ) 





co 


Generate( float(range =-1 ..1, digits = 10) ) 





co 


Generatel float(range =-1 ..1, digits = 10) ) 








co 


Bx := BE 1] : By := BI 2] : B := 13] : 

















GS := GramSchmidt({œl, «2, «3}, normalized) ; 
1 0 0 0 


0  GS[1][1] GS[2]1[1] GS(3]1[1] 








“on 0 GSI1]1I2] GSI2]1I21] GSt3]12] 
0 GSt1]1I3] GSI21131 GSt3113] 

M :=A:Q: 
B 
return M, M%T-G°M, À QQ-QT, | Py 
B 


end proc: 





1 
= : 
Y Bx°7 By: Be 
Bx-y + 2 PE RAR 
A := 
By . BB 14 — Be _ Tr B& 
Br D RR ER 14e 
Generate(float(range =-10 ..10, digits= 10) ) 

ol := | Generate(float(range =-10..10, digits=10)) |; 
Generate(float(range =-10..10, digits = 10) ) 
Generate(float(range =- 10.10, digits = 10) ) 

2 :— | Generate(float(range =-10..10, digits=10)) |; 
Generate(float(range =-10..10, digits=10) ) 
Generate(float(range =-10 ..10, digits = 10) ) 

œ ‘= Generate(float(range =-10..10, digits=10)) |; 
Generate(float(range =-10..10, digits=10) ) 





pl(); 
1.37137909400000  0.0693759093241285 —0.867398928126097  —0.351406747366782 
0.407386895800000  0.749462011891991  —0.292339217513503  —0.720283663541724 
0.617131505500000  0.240561552388650  —1.09680276211141  0.346417572003398 
0.577810731800000  —0.620684982570698  —0.681133160210365  —0.696184609664106 


1.371379094  0.4073868958  0.6171315055  0.5778107318 
0.4073868958  1.069986314  0.1060190203  0.09926397720 
0.6171315055  0.1060190203  1.160603294  0.1503703932 
0.5778107318  0.09926397720 0.1503703932  1.140789485 


#: 


1 0 0 0 
0 0.7375437000  —0.143325946099933  —0.659914360441726 
O0 0.2225070544 —0.871069471880707  0.437868228763262 | 
0 —0.6375891314 —0.469782553941856  —0.610560767365580 
1. 0. 0 0. 
0. 0.999999999352344 2.12288964185348 10-10  2.59283428061252 10-10 REVERS 
0.4500079578 
0. 2.12288964185348 10-10 0.999999999860575 —_1.43071110514370 10-10 
0.4213355259 
0. 2.59283428061252 10-10  —1.43071110514370 10-10 0.999999999113585 





#(2)Maintenant un programme quigénére 2 matrices de Lorentz symétriques et leur produit qui n'est pas symétrique : 


p2 :=prod() 
local M1, M2, M ; 


M1 := pi( 131; 


M2 := pI( )131; 
M := MI «M2; 
return M1, M2, M; 
end proc: 


pA(); 

1.368682988  0.3848997155 0.7760112340  0.3506449710 
0.3848997155  1.062544372  0.1260981333  0.05697813945 
0.7760112340 0.1260981333 1254231334  0.1148758352 
0.3506449710 0.05697813945 0.1148758352  1.051907282 


, 


1.705641095 —0.8360995954  0.8859238484 0.6521410478 
—0.8360995954 1.258372234 —0.2737689683  —0.2015252013 
0.8859238484 —0.2737689683 1.290083214 0.2135343479 | 


0.6521410478 —0.2015252013  0.2135343479 1.157185647 


0.999999999726752 1.35069122553233 10—10 —2.73939648742783 10-10  —7.82773845742213 10-12 
1.35069122553233 10-10 —0.999999998551066 1.71263558890189 10-10 —2.41488051777594 10-10 
—2.73939648742783 10-10 1.71263558890189 10-10 —0.999999998495485 —3.83972298401147 10-11 
—7.82773845742213 10-12 —2.41488051777594 10-10  —3.83972298401147 10-11 —0.999999998796079 





(D) 


(2) 


2.92882429158172 —0.943119770924961  2.18316911404290  1.38647374570310 
—0.0830210202213288  0.969257551838184  0.224944035647489  0.129740003136537 
2.40423475039461 —0.856664283201202  2.29555782999394  0.881410965250509 
1.33819823203237 —0.464909675392575  0.667863617876428  1.45896939303546 


#(3)On considère la HERO polaire : M = A: Q qui correspond aux changements de base 
#M:808" & MB Bo 38" X=M-X'. 

#Ce qui revient à réorienter les axes Ox'y'z' selon les axes Oxyz 

#puis d'appliquerl'opérateur symétrique ad hoc dans cette situation . 


Un programme qui donne les parties symétriques et orthogonale d'une matrice de Lorentz 
#par la formule :Q = A"1M. 
# En entrée on a la matricede Lorentz considérée . 


#En retour on a B, les parties symétriques et orthogonale de la matrice la vérification que 
#Q est bien orthogonale . 


p3 :=proc(M) 
local B, y Bx, By, Be, Q, Aoinv; 
7:= MU] 1] : 
MT21111 M131111 MT41111 
BEI] = :B21 = ——  —— :B131 :— 


7 ‘à 7 
Bx := BL 1] : By := BI 21 : Be := BL3] : 








Y -Bx-y —Bp-Y By 
À À 
HONTE +1 +1 “Bc-fy +1 Bc- ke 
Apinv := Ê Ê Ê 
hr Sir bR IDR Sr BR 
ê À Ê 
_ —— —— .p. —— .p2 
DR Re LR ASE 
Q := Aoinv - M : 
Bx 
return | By , Aginv, Q, Q + OT; 
Be 
end proc : 
M = p2()t31; 
P3(M); 
1.38083390331948 0.932012422794345 —0.143386261535510  —0.132270532031612 
0.891753407171495 1.32442538856998 —0.125725949457531 0.159105358690483 
M := 
—0.0627733200163432  0.0563998342199275 0.999893227938980 —0.0311942286970146 
—0.327929322210753  —0.333711858912975  0.0704708246877686 0.995594244327785 
1.38083390331948 —0.891753407171495 0.0627733200163432 0.327929322210753 
0.645807873798396 
—0.891753407171495 1.33401075904255 —0.0235120652163062  —0.122827589940295 
—0.0454604423207151 | k (3) 
0.0627733200163432  —0.0235120652163062 1.00165508803465 0.00864621940118425 
—0.237486435857652 


0.327929322210753 —0.122827589940295  0.00864621940118425 1.04516805654341 








0.999999999283068 4.29898894260816 10-11 —4.47618851706633 10-11 —3.12159520454713 10-10 
2.68530746250217 10-10 0.975335411132845 —0.0720198983434591 0.208647957072290 
—1.89027157641852 10-11 0.0809733729268503 0.996112690509814 —0.0346816876264137 
—9.87482318137722 10-11 —0.205339106119544 0.0507212069135954 0.977375675632391 

0.999999998566135 2.48552607731753 10-10 —4.91833422936703 10-11 —4.14943236258999 10-10 
2.48552607731753 10-10 0.999999999957519 —5.79414641266318 10-12  3.58315876969328 10-10 
—4,91833422936703 10-11 —5.79414641266318 10-12 0.999999998774467 —7.61509605373689 10-12 
—4.14943236258999 10-10 3.58315876969328 10-10 —7.61509605373689 10-12 1.00000000065062 

#(4)Un programme qui donne les parties symétriques et orthogonale d'une matrice de Lorentz 
1 0 0 MT2]1(21] MI2]1[13] MI2][4] 
Y 
#par la formule :C-1.71= 0 1 0 |- FPS LS MI31121 MI3113] MI3][4] 
Y 
0 0 1 MI41121 MI41131] MI4]114] 


# En entrée on a la matricede Lorentz considérée . 


#En retour on a la partie orthogonale Q de la matrice M, la vérification que 


#Q est bien orthogonale etque M = A:Q. 
p4 :=proc(M) 
local B, Bx, By, B, y Cinv, N, A, Q: 
unprotect(}) : 
Y:= MTL, 1]: 


M2, 1] 








MT2]1121 MI2]113] MI2114] 


1 0 0 
Y 
N:= 0 1 0 | —.p.prr|| MI31121 MI31131 M(31(4] 
0 0 1 MT4]1121 MT41131 M141141 


1 0 0 0 
0 NL1I] NI12] Mi1,3] 
0 M21] M22] M2,3] 
0 M31] M32] M3,3] 
Be := BIT : By = BL21 : Be := BL3] : 














Y Bx-7 By-7 B°Y 
ê ê À 
: ——  py2  ——.By. — — pr. 
Bey 14 Be rh or BR 
A := : 
Ê À À 
BE er DU PSS or D 
À ê À 
| —— .p. —— .p. —— .p72 
By Gr Rk& er kh 1455 & 
return "] =" M, B Q GOT, À - 0: 
end proc: 
P4(M); 
1.62544210582349 —0.620767627885989  —0.531312165942930  —0.987125611695587 
"M=" —0.742738427361049 1.16928792128352 0.202238837494802 0.378847701738497 
? | —0.498020524205100  0.0655819952579160 1.09942281087722 0.187063969441193 | 
—0.917811105247344 0.117546614876760 0.180725079895222 1.34011136198897 
1 0 0 0 
—0.456945482524434 
0 0.993672557106914 0.0519304525664274 0.0995895445522121 
—0.306390810488319 
0 —0.0521715132503900 0.998638122894378 —0.000184042385260623 
—0.564653211553392 
0 —0.0994634731377091  —0.00501285935373779 0.995028586607148 
1. 0. 0. 0. 
0. 1.00000000003526 —4.98607929168846 10-12  1.02634983845107 10-11 


0. —4.98607929168846 10-12 


1.00000000116444 


1.69018843788459 10-11 


, 


0. 1.02634983845107 10-11 1.69018843788459 10-11 0.999999999412934 
1.62544210582349 —0.620767627787629  —0.531312166142487  —0.987125611643738 
—0.742738427361049 1.16928792131431 0.202238837521153 0.378847701787455 
—0.498020524205100  0.0655819952785598 1.09942281089489 0.187063969474021 
—0.917811105247344  0.117546614914805 0.180725079927785 1.34011136204947 





t 
#(5)Transformation d'une matrice de Lorentz M symétrique sous la forme M = A :* Ao.A : 
#D'abord un programme qui construit une base orthormée à partir de B : 


p5 :=proc( P) 
local A, ci, c2, c3, di, d2, d3, c, d; 
with(Random Tools) : with(LinearAlgebra) : 


randomize ) : 
if Norm(B, 2) =0 then print("B=0");return send if: 


cl = Gencrate| float(range =-1 1 )); 
e2 = Gencrate( float(range =-1..1 )); 


e3 = Gencrate| float(range =-1..1 )); 


d1 := Generatel float(range = 


-1.1)); 


da2 := Generate( float(range =-1 1 )); 


d3 := Generate| float(range =-1..1 1); 


cl d1 
le ler: 
c3 d3 


À := GramSchmidt(| 6, & d], normalized) : 


1 0 0 0 
0 AI] AZ] ASIA 
Flo A2] A2] ABJI2] 
o AU]I31 4213] AS]13] 
end proc: 
#le programme final : 


#On donne la matrice de Lorentz M symétrique en entrée eton obtient À, Ao et 


t 
#A + Ao + A=M, la décomposition n'est pas unique. 








p6 :=proc(M) 
local} B, À, Ao : 
unprotect(?) : 
#calcul de A 
7:= MI1111; 
MT2111] 
Ÿ 
M131111 
B=| —— |; 
Ÿ 
MT41111 
Ÿ 
A = pS(B) : 
#calculde Ao 
Ÿ y-Norm(B2) 0 0 
je y Norm(B, 2) 7 0 0 ’ 
0 0 1 0 
0 0 0 1 


return À, Ao, À * Ao.AT : 
end proc: 


#programme pour transformer une matrice de Lorentz M nonsymétrique sous la forme M = À + Ao -B, 
#il donne À ,/0,B ,A + Ao -B : 


p7 :=proc(M) 


local}, B A, B, Ao, x, By, Br, © N': 
unprotect(?) : 

#calcul de À ,Ao 

À := p6(M)[ 1] : 0 := p6(M)[2] : 
#calculde Q : 
Q := p3(M)[3]; 


#calcul de B: 

B := AT. Q : 
N':=A:°/A0:-B: 
return À, A0, B, N': 
end proc : 


#A:=p2( )L11;p6(A); 











M 2= pl()tn; 
P7(M); 
1.14630177500000 —0.187578579359841 0.456368299099120 0.265612520527367 
0.183822213000000 0.394658776379304 0.851535493083302 —0.391052887222597 
— —0.463019536400000 0.104666235055864 —0.576945157041309  —0.933041454534980 
—0.256573697000000 0.931920142736560 —0.387677203626374 0.216936152945549 
1 0 0 0 1.14630177500000  0.560363952331461 0 0 
0 0.328040753250266 0.0775274998606912 0.941476898823284 0.560363952331461  1.14630177500000 0 0 
0 —0.826283586795056 —0.459507533735608  0.325742629317178 | 0 0 110} (5) 
0  —0.457869739758195 0.884783765793880 0.0866774968671706 0 0 0 1 
1.00000000029561 —1.63093039073914 10-10  5.90026638658259 10-12 —7.45295700155246 10-11 
—7.71801484377995 10-11 —0.334744193368859 0.814414091589783 0.474000012971911 
—6.14795999036710 10-18 0.807049798004479 —0.0118824321028395 0.590363811326780 j 
6.17397736510085 10-17 0.486432880739307 0.580162473236776 —0.653295153483057 
1.14630177529561 —0.187578579403134 0.456368299204451 0.265612520588671 
0.183822213042427 0.394658776292668 0.851535493202813 —0.391052887232183 
—0.463019536506866 0.104666235532494 —0.576945157289849  —0.933041454418524 
—0.256573697059218 0.931920142085773 —0.387677203752586 0.216936152344186 
1.241716459 0.3868137862 —0.01793713486 0.6260296461 
0.3868137862 1.066745687 —0.003095097519 0.1080229825 
— —0.01793713486  —0.003095097519 1.000143524 —0.005009187556 (6) 
0.6260296461 0.1080229825 —0.005009187556 1.174827247 
FJ 0 0 0 1.241716459  0.7361112442 0 0 
0 0.5254827843 0.214368773706647 —0.823355252697488 0.7361112442  1.241716459 0 0 
O0 —0.02436742407  —0.963548601974385  —0.266421320706501 | 0 0 110} (6) 
0 0.8504552145 — 0.160062863999307 0.501104587495512 0 0 0 1 
1.24171645900000 0.386813786156753 —0.0179371348501167 0.626029646081973 
0.386813786156753 1.06674568698906 —0.00309509774445957 0.108022982466224 
—0.0179371348501167  —0.00309509774445957 1.00014352415918 —0.00500918760113681 
0.626029646081973 0.108022982466224 —0.00500918760113681 1.17482724744993 
1 0 0 [1 0 0 
Ê Y 
#(6)vérification de la formule : 0 1 0 |+ Fer -BeBT| =|0 1 0 \- Fe -B+B#T; 
0 0 1 0 0 1 
unprotect( y) : Digit= 16 : 
1 
B:=pi()té:y:= ; 
1 — (Norm(B, 2))? 
—0.008610915414 
be —0.3926553995 
—0.4894613276 


7:= 1.284401838 


(7) 


-1 





1 0 0 1 0 0 
0 1 0 P 0 1 0 di 
— —— .B.BT| :-h:=—  ——— .8. BST. 
0 0 1 0 0 1 
0.999958310463802 —0.00190103150520501  —0.00236971503647201 
a:= | —0.00190103150520501 0.913313480710222 —0.108058360767806 
—0.00236971503647200 —0.108058360767806 0.865300747711506 
0.999958310463786 —0.00190103150594011  —0.00236971503738834 
b:=| —0.00190103150594011 0.913313480676701 —0.108058360809590 
—0.00236971503738834 —0.108058360809590 0.865300747659420 
1 0 0 
re . 1 0 Ê 0 
#(7) Vérification de + Toy -B-B*T| -B=yB: 
0 0 1 
1 0 0 
0 1 0 + ge gr By BR 
1+7 U 
0 0 1 
—0.0110598755850366 
—0.504327316838143 
— 0.628665028823940 
—0.0110598755846041 
—0.504327316818424 
—0.628665028799360 








#(8) Vérification de la formule (cas du produit de 2 matrices sans rotation ) 


y B°B%T , 
#0 | —— | (y-72-p1 « B2%T + CI - C2) : 
1+7 
unprotect(?) : 
1 0 0 
Id:=| 0 1 0 
0 0 1 


a = pl( ) : MI := al3]; BI := al 6] : 1 := 


; YL.BI =, À - BI; 


1 — (Norm( BI, 2))? 
1 





5 V2.B2=,2-B2; 
1 — (Norm(P2, 2))? 


a = pl) :M2 := al3]; f2 := al 6] :2 = 


M := MI *: M2; 


= A 2(BPA RE +1); fe MU, D 





= . . T « 
= y BI - BIST: 





= . ° B2%T : 
QG + BB: 





(8) 


(9) 




















M2, 1 
MI, 1 
"1 - BI + C1 -B2 M3, 1 
BE per g+nn | un L 
M4, 1 
M1, 1 
"y-B = 4B 
M2, 1] 
"y : BI — M3, 1] 
MT4, 11 
a 
C:= |14+ 1+7 
v-B:B#T 
#—= _ à . + B2%T . . 
a (ua Iéy (91-72 * BI + B2#T + CI * C2); 
QI := p3(M)13]; 
"CQ =", ca; 


"vérification de l'orthogonalité de Q", QT + Q; 
Bloc := (-2+B1-B2#7+C1-C2); "M=", M; 


1.094952487 —0.2539834648 —0.2300330515 0.2854787980 
—0.2539834648 1.030791916 0.02788826564 —0.03461028101 
MI := 
—0.2300330515  0.02788826564 1.025258427 —0.03134656248 
0.2854787980 —0.03461028101  —0.03134656248 1.038902144 
YL := 1.094952487 
—0.253983464845897 
pi = | —0230033051513705 
0.285478797969767 
1.638170280 0.8327160307 —0.5379022638  0.8371660729 
0.8327160307 1.262839739 —0.1697842786  0.2642443568 
M2 := 
—0.5379022638  —0.1697842786 1.109674060 —0.1706916074 
0.8371660729 0.2642443568 —0.1706916074 1.265656481 
Z := 1.638170280 
0.832716030705982 
72.B2 = —0.537902263803150 
0.837166072928481 
1.94495098291234 0.705536233442017 —0.849845567264778 1.25012617845242 
0.398313074866063 1.07634833098581 —0.00153941170382430  0.0111895876706808 
M := 
—0.931341409991784  —0.338689714120093 1.26205358646206 —0.399883538410058 
1.32543736156987 0.473861715598860 —0.359799754262568 1.55009141958829 


7 := 1.944950983 
Ÿ := 1.94495098291234 


0.204793374191607 
B:=| —0.478850838892369 
0.681475971905856 





0.204793374416889 
—0.478850839005314 
0.681475971998628 


Il 


0.398313074445852 
"Y*B =". | —0.931341409814088 
1.32543736144917 


0.398313074866063 
"yL* BI =" | —0.931341409991784 
1.32543736156987 

0.980922529528066  0.113404640943630 —0.157893566680858 
Q:=| —0.115563721858629  0.993289735023370  —0.00453085372986185 
0.156320239142954  0.0226911845224586  0.987445741923861 
0.999999998911629 —5.99240324028472 10-10  1.17939602528594 10-10  6.07358607851438 10-10 


1.47205275746742 10-10 0.980922530025584 


0.113404641007299 —0.157893566935660 
Q1 : 


Ï 


—3.44197559343229 10-10 —0.115563721969872 0.993289734916023 —0.00453085341501858 





4.89844609319334 10-10 0.156320239228924 0.0226911845669175 0.987445740927857 





1.07634833046138 —0.00153941164370804 


0.0111895876454864 | 
"CA =",| —0.338689714138731 1.26205358651223 


—0.399883538411865 
0.473861715729400 —0.359799754263803 1.55009142019491 


"vérification de l'orthogonalité de Q"! 


0.999999999911268 —1.76368924861377 10-11  1.36538543895174 10-11 


—1.76368924861377 10-11 1.00000000014538 —2.01605528363526 10-10 


1.36538543895174 10-11 —2.01605528363526 10-10 1.00000000027829 


1.07634833045804  —0.00153941163968860 
Bloc := | —0.338689714130929  1.26205358650283 
0.473861715718296  —0.359799754250428 


0.0111895876395738 
—0.399883538398040 
1.55009142017523 
1.94495098291234  0.705536233442017  —0.849845567264778  1.25012617845242 
nn | 50706606 107685500821 


, 


—0.00153941170382430  0.0111895876706808 10 
1.26205358646206 —0.399883538410058 (D 
—0.359799754262568 1.55009141958829 


—0.931341409991784  —0.338689714120093 
1.32543736156987 0.473861715598860 


Annexe : 


Démonstration du Theorème 3 : (J-M. Souriau."Calcul Linéaire " - PUF 1964 p.378). 


Soit X € R' tel que XX =1 Soit N la matrice définie par : 
©O : étant la matrice nulle de M CR) ; 
h — h — 


0 X 
N= x © . Ona alors Va ER: 
n —1 
ch( @) sh (æ) x 
M = exp\( aN) = 
FN lt (4, _, +(eh( a) -1)XX) 
Ri — 1 
est une matrice de M (R) deLorentz de déterminant égal à 1 . 
n 
Démonstration : 
1 0 
Ona N° = Ft t en notant 0" — À la colonne nulle d'ordre n — 1, 
0 XX 


puisque IXX=1et ona N°=N. 
En remarquant que pour toute matrice carrée À et tout & € R, 

















C4" + œ RE + œ de 
ep(A)= Der » sh(a) = 2 Cn+1)! : ch(a) = 2 (2m)! 
. se 
a a a 
donc exp(aN) GE SEINE 21 + 31 Æ 11 RÉ Gusss + 
2p—1 2p 2 
+ a N LUN Fi 
2p—1)! DT CPE 
=, +sh(œ)N + (ch( @) -1 N°. 
Vérifions que M est bien de Lorentz : MGM = G ? 
On remarque que M est symétrique : 
ch( @) sh (@)'X 1 ‘0 
L ; 0 - Id 
sh (@)X (14 , + (ch @) 1)XX) ne 
ch( @) < sh(æœ)'X 
— t 
sh(@)X (1, _,+(ch(a) -1)XX) . 
ch( @) - sh(œ)'°X 


a b 
c d 


a = ch” a) — sh (æœ) ‘XX =] ,car XX =1. 


.Calculons a, b, c, d. 








b=ch( @) sh(@œ)'X- sh(a@ a)'X (I, _, +(ch( a) - 1)XX) 
= ch( @) sh(@œ)'X- sh(œ)'X- sh(œ)(ch( æ) -1)'XXX 
= ch( @) sh(@œ)'X- sh(œ }x X - sh(aæ )ch( æ)'X X +sh (æœ )'X =0, 


c=ch( @) sh(æœ) X- (4, _, + (ch( a) L)XX }sh (a) X 
=ch( @) sh(@œ)X- sh(œ)X- sh(æ) (ch( æ&) -1)X=0, 


+ (ch( æ) -1 XX) (4, + (ch( æ) -1 XX) 


=, (cha) -1)XX + (ch à) -1 XX + (ch( @) -1 VXXXX, 


on remarque comme XX =1 :X CXX)X =X'X (association du produit matriciel) , 


et(ch( æ@) -1) =ch?( æ) —2ch( @) +1 
“HW, + (2(ch( @) -1) +ch°( @) —2ch( æ) +1)XX 


= Id + (ch ( @) —1)X'X = Id +sh( œ)X'X 
m1 — 1 1 — 1 























finalement : 
2 => 
h° (a)XX (4, + (ch( @) - xx) D R. 
b 1 0 
Donc — 
d 0 EC 
et M est bien de Lorentz. 
Une autre démonstration possible est de calculer : 
ch( @) sh(œ)'X || 7 
aN)G= 
FPieNeS | (a, _,+(eh( a) -1)XX) | 0 
R -1 
ch( æ) -sh (a) x 


et 
— 1 








sh (@œ)X - (14, + (ch( @) - 1)XX) 





0 hs -sh (œ)X (14 , + (ch @) -1)XX) 
ch( æ) -sh(æ@)'X 
| sh(a@)X (Hd, +(ch( a) -1)XX) 


donc exp( aN)G= Gexp( -aN) —  exp(aN) Gexp( aN) =-G 
comme M=exp( aN) est symétrique MGM = G et M est de Lorentz. 
Pour une troisième démonstration voir : J.M. Souriau : Calcul lineaire.PUF. 


(4) 


Montrons rapidement que det(exp(A)) =" pour A EM, (R) : 


Considérons À comme élément de M, (©). À est alors trigonalisable et peut s'écrire : 

A=P BP avec B trian gulaire supérieur dont la diagonale est composée 

des valeurs propres À, de A. 

en s'appuyant sur la définition de l'exponentielle on peut écrire : exp(A) =P lexp (B)P. 

En développant exp(B) en série et en remarquant que B° triangulaire supérieur V k > 1. 
À 


Alors exp(B) a pour diagonale les e À 


22 


+ 
i 


À 
D'où det(exp(A)) =det(exp(B)) = [ le =ef =e"t4) 
i 
Comme aN est de trace nulle det( exp( aN) ) =1. 


Lemme 2 : 
Si M est une matrice de Lorentz dont la première colonne K, est 


de la forme: K,= | a, 0,.., 0) alors M est de la forme : 


t 
€ 


M= où e=+1,"0=(0, 0,0), QQ=Id (LD). 
Q a ou 
Réciproquement toute matrice de cette forme (1) est de Lorentz. 
Démonstration : 
. œ ‘L t 
Soit M= avec CEM(R), L={(l,1,1) 3 
0 C ; 
Comme G)M= G on a: 


M 

æœ 0 a ‘!L 
G= 
L 


C 0 C 


2 
a a'L 


oL L'L—'CC 


par identification œ=+ 1, L=0, "CC cl. _ 


Si M est de la forme (1) un calcul immédiat montre que MGM =G . 








Theoreme 2 : 
Toute matrice M de Lorentz peut se mettre sous la forme : 
£ 0 ; 0 x 
M = exp( aN) ouaER,N= , À 
0 @Q X O 
e R’' — 1 
telque: 'XX=1 , ‘OQ=Id,_, 
Démonstration : 


Soit Mune matrice de Lorentz M= (M ;) et M, la première colonne de M. 

Si M, = (M; 7 ) alors : 
t _ 2 2 2 =. 
M,GM, = Mig My ee Min © G,1=1. 

On peut poser M,= (B Y) avecB=m, ,et = (M 290 Min). 

On a donc — {YY=1 donc|Bl > 1.On peut poser B=£-ch(æ) avec 

£=tlet a ER. 

fs =. 
Deplus YY= f —-1=5h (æ). 


Si à # 0 posons jee, on a ‘XX=1I. 
£-sh( a) 


Siæ=0 ona M,=(8& 0,0,0). 

Danns tous les cas M, =£(ch(æ), sh(æœ)X) avec XX 1. 
d'après le lemme 1 M, est la première colonne de £-exp( aN) avec 
0  %X 
X O 


On remarque que exp( - aN) -M est une matrice de Lorentz puique M et 
exp( - aN) le sont .Calculons la première colonne K , de ce produit : 


N= 








ch( @) (sh(æ))'X 


KT | (sh (a@))x 1, +(eh( a) -1)XX 


- M, 


En remarquant que ‘XX =1, ch°(@) F sh (æ) =1 ,etque 
(4, + (ch( æ@) -1)XX ) -X=ch( @)-X 


on obtient : K, =(£,0,.., 0 ).On peut alors appliquer le lemme 2 et 
le théoreme 2 est démontre. 
Théorème : 
La décomposition du théorème précédent est unique . 
Démonstration : 
Elle découle des resultats qui suivent : 


Lemme : 
Toute matrice carrée symétrique À d'ordre n > I, à coefficients réels , est diagonalisable . 
De plus les sous — espaces propres de À sont orthogonaux 2 a 2. 


On peut donc écrire que À = PD fP avec P matrice orthogonale , D matrice diagonale. 


Démonstration : (D'aprés J. Grifone Algèbre Lineaire .Cépadues éditions 2002) 


: a. C4 — 
On note, si A -(hoasne alors A rer 


et si À est symétrique ‘A =A. 
(a)Montrons que toute valeur propre de À est réelle. 


On peut considéèrer À comme un opérateur de ©' dans ©! et 
P ,(X)le polynôme caractéristique de À .Comme 


polynôme sur ©, il admet au moins une racine À, montrons 


que cette racine est réelle Soit IX = (x sXoputes Xy) Æ0, 


le vecteur propre associe à À. 

Comme À est symétrique on a , en notant À et X les conjugués de À et X dans C: 
((AX)X = XAX puisque À est symétrique .D'autre part AX = X d'ou : 
comme À est réelle AX = ÀX, il vient '( AX)X = IXAX c'est à dire que : 


AD xl =22 xl? doù i=1. 
i=1 i=1 


Comme AX=AX = AX=AX = AX=AX = A(X 4X) fx +x) 

on peut donc associer un vecteur propre réel à À .Ou plus simplement remarquer que 
À est aussi valeur propre de À comme opérateur de R' dans lui — même. 

Cela implique que tout endomorphisme f, d'un espace vectoriel réel E de dimension n 
dans lui — même, symétrique , admet au moins une valeur propre réelle À 


etun vecteur propre réelx associé : f(x) = x. 
Rappel : f'est dite symétrique si(f(x),y)= (x, f(y)) VxEëE ,VyEE <,> étant 
le produit scalaire euclidien de E. 
(b)Montrons par récurrence sur n, dimension dE, 

qu'il existe une base de n vecteurs propres réels orthogonaux 2 a 2, 
pour tout endomorphisme symétrique f d'un espace vectoriel réel E 
de dimension n dans lui — même . 

Sin=liln'ya rien à démontrer puique dans ce cas f=À. 
Supposons que la propriété soit vraie au rang n — I : 
Pour tout endomorphisme symétrique f d'un espace vectoriel réel E 
de dimension n — 1 dans lui — même, il existe une base de n — 1 
vecteurs propres réels, orthogonaux 2 a 2. 
Montrons qu'elle est vraie au rang n. 


On vient de montrer qu'il existe À etx valeur propre réelle 

et vecteur propre réel associe pour f endomorphisme symétrique 

d'un espace vectoriel réel E de dimension n , n quelconque . 

Considérons le couple (x, À) associé à E de dimension n. 

Considérons dans E, H=(x { ) l'hyperplan orthogonal de x dans E , 
de dimension n — I. 

Montrons que f'est stable sur H : 
Soit y un vecteur de H : (f(y),x) = (y, f(x)) =À{x, y) =0 donc f(y) est orthogonal a x. 
Soit y etz 2 vecteurs de H, comme ils sont aussi vecteurs de E on a aussi (f(x), y) =(4f(y)) 
et la restriction de f à H est aussi symétrique . 

Comme H est de dimension n — 1, l'hypothèse de récurrence s'applique : il existe une base 
de vecteurs propres, dans H, pour la restriction de f a H : A ee , X Orthogonaux 2 2: 


AÏOTS X,X3s us. , X, est une base de n vecteurs propres réels de f dans E : 
f'est diagonalisable .Comme H est orthogonal à x, alors x,, .…… , X, sont orthogonaux à x. 


LES X,X,, se ,X, , apres normalisation , forment une base orthonormée de 


n 
vecteurs propres Soit 8 une base de E , À la matrice qui représente f dans 8, 

P la matrice composée de colonnes représentant successivement les composantes des 
X 3 Xe, X, , D la matrice diagonale composée des valeurs propres À et À, associées à 


X y Xp es X,. On a donc À = PDP ! de plus P est une matrice orthogonale 


et done P°|="P. 
(Pour d'autres démonstrations consulter J.M. Monier Algebre 2 Dunod 1997 par exemple) . 
Remarque : La propriete est fausse si les coefficients de la matrice sont des nombres 


2 t 


complexes .Considérons: A=| ñ 
i 
Son polynôme caractéristique est X?-2X+1- (X — 1 a 
Si la matrice À était diagonalisable , elle serait semblable à la matrice unité. 


Lemme 
Soit E un espace vectoriel de dimension n sur R 
Soit fet g 2 endomorphismes symétriques de E qui commutent :fog=gof 
alors il existe une même base de vecteurs propres dans laquelle f et g se diagonalisent . 
Démonstration : 
Puisque f et g sont symétriques ils se diagonalise chacun dans une base orthonormee de 


les espaces propres correspondants Montrons que E est stable par g : 
i 


Soit x CE; : f(x) =Àx alors f(g(x)) =g(f(x)) =Àg(x), 


ou bien g(x) =0 ou bieng(x) Æ 0 etg(x) est valeur propre f 
et dans les 2 casg(x) € E, donc 8(Fa) CE, . 
i i i 


On peut definirg,yx :E 457 E n tel que g Æ, soit un endomorphisme symétrique 
L L 
i i 
donc il existe une base B ortonormée de E a constituée de vecteurs propres de f 
L 5 
L 


qui diagonalisent gr 
i 
P P 
Comme D E, =E, U B = base de E composée de vecteurs propres de f et g qui 
1 i=1 ! 


i=1 i 


diagonalisent f et g . 


Lemme : 

Soit M une matrice inversible à coefficients réels alors il existe un couple unique de matrice 
O orthogonale et S symétrique définie positive telles que M = OS. 

Démonstration : 

Existence : 


MM est une matrice symetrique définie positive car { UM ) = MM, 

det( MM) = (det(M) + > O0 car Æ 0, Métantune matrice inversible, 

de plus X'MMX = MX || : > 0 avec || - || norme euclidienne , 

UM étant symetrique, elle est diagonalisable Si À est une valeur propre de MM, 

elle est différente de 0, M étantune matrice inversible et soit X a 

un vecteur propre associe on et si À < 0 on aurait || MX || 2 =X, MMX, = AX,X, < 0, 


ce qui est impossible, donc on peut écrire : 


À, 0 
t _ N -1 th nl 
MM=O0, . O, avecÀ, > 0 car O=0O . 
0 À, 
JA, 0 
. N * . fe. 
SoitS=0, | O, ,Sestdelaforme S=0, DO, , 


ue - : 
S'est symétrique : s=0 D 0, =0, D 0; =S et S'est définie positive . 
Soit O=MS | .Montrons que O est orthogonale : 


-] 
On a O=(s!)m=('s) M=S ! M Donc : 


‘00=S ! (MM) S'! = 


O0 è 0;'0 . 
1 \ 1 71 \ 
(a) | L 
= 
(72) 0 
0j'0) a 0j! 
-1 
0 (/2,) 

= Id. 

donc O etS répondent à la question . 

Unicite : 


Soit P un polynôme tel que P(1;) = JA, pouri= 1, p 


par exemple le polynôme d'interpolation de Lagrange pour la fonction sur R+ :f(x) =Vx , 








{x —2 
car [I j) 1six=let-0six=21 #1. 
: Li ve i j i 
RATE 
À, 0 
t à -1 
P(MM) =P\0, . 0, 
0 À, 
À, 0 
N -1 
= O,P . 0,'=-S 
0 À 


SiM=OS alors M=S0 ‘et MM=S° 

soit S,une autre matrice symétrique , définie positive et O,une autre matrice orthogonale 
telle que M=OS, alors comme précédemment _S ra =MM . 

Calculons S, MMS=S, se —S; 2, = MMS, donc S,commute avec MM et 


commute donc avec P('M M ) et donc avec S .DoncS et S Sont diagonalisables 
dans une même base (lemme 4).Leurs valeurs propres sont positive non nulles de même 
carré carS,? ='UM=S* donc elles sont égales etS, =S et donc O = O,. 


Corollaire : 

Soit M une matrice inversible à coefficients réels alors il existe un couple unique de matrice 
O orthogonale et S symétrique définie positive telles que M=S 0. 

Démonstration : 


On considère N=M | = ON Sy de manière unique d'aprés le lemme3, 
_. ete \ et eee 
S,, etant symétrique definie positive, \Sy ]=\Sy) =Sy 


donc Sy ‘est aussi symétrique donc diagonalisable . 
Soitx Æ O un vecteur de Se et À Æ 0 la valeur propre associée : 


-1 -1 
ona X=Sy-Sy -x=Sy-À-x etdoncS, x est vecteur proprede S,, , 


etu= nr est valeur propre de Sy, donc À > 0 et Se est symétrique définie positive. 


Comme O,, étant orthogonale Oy = 0” : ( est aussi orthogonale . 
On peut donc écrire M =$, 0)! de maniere unique. 


On peut aussi considérer N ='"M = OYSYy, et donc M='SY0y = SyOy. 


Cas où M n'est pas inversible : 
On remarque tout d'abord que GL, (R) est dense dans M, (R) : 


Soit M E M, (R), A, = min ( A; / À, valeur propre de M ) et considérons les matrices inversibles 
M — Al, , avec0 < À < À,, qui tendent vers M lorsque À — 0. 


On peut donc considérer une suite( M;) de matrice inversibles qui tend vers M. 
D'aprés ce qui précède M, = OS; avec O; orthonormale et S; symétrique . 


On remarque O,(R) = {O | O orthormale }est compact dans M, (R) : 
1 1 


Si0 EO,(R) et p le rayon spectral IO1,=p° (0/0) =p" (jy) = /n, 
donc est borné dans M, (R) D'autre part soit D: M — MM définie sur M, (R), 


-1 
qui est une application continue dans M,,(R) et on remarque que ® (4) =0,(R) 


n 
par definition même d'une matrice orthormale .Donc O,(R) est ferme dans M,(R) 


et donc compact. 
De la suite des (0;) on peut donc extraire une sous — suite convergente (Si ) convergente 


vers la matrice O qui est orthonormale car O,,(R) est ferme dans M, (R). 


On en déduit que la suite de terme général; = ‘0, M. est convergente 
k k k 


vers une matrice S symétrique (considérer les composantes) et 

M=OS avec les proprietes requises. 

Puisque O est inversible le rang de S est égal au rang de M, que l'on suppose inférieur 
strictement à n.O n'est detemine par M que sur Im(S) d'ou l'impossibilité de démontrer l'unicite. 


Remarque relative aux droites d'univers : Rappel de géométrie affine : 

Lemme : 

Toute intersection de sous — espaces affines est un espace affine, éventuellement 
vide et, sinon, de direction l'intersection des directions des sous — espaces 

des sous — espaces considérés . 

Conséquence : 

L'intersection d'une droite affine D et d'un hyperplan affine d€ est un sous — espace 
affine de dimension au plus égale a I : 


Soit D est parallele à 4 soit D() H=9 ou DC. 


Si Dn'est pas parallèle à 4€ D N Ge se réduit à un point unique P. 


Définitions : 
On considere l'espace vectoriel réel E de dimension 4 de la relativité restreinte. 
On considere les hyperplans affines æ ={(a,x, »,27)/(x, 7,27) ER}. 


Ils représentent l'espace physique de l'observateur à un instant t donné ( présent absolu). 
On considere les droites affines de l'espace E qu'on va partager en 2 categories : 


= {A droite affine/21 # telque À C #} et 





dy 
C X) 
C;= D> droite affine d'equation parametrique :t| _, | + 
"70 4 Yo 
9) 
y 
>. X9 
V etant un vecteur de dimension 3 etX, = un vecteur constant. 
Vo 
79 


A est une droite d'univers dont tout les points sont contemporains , D> est une droite d'univers 


qui représente la trajectoire d'un mobile animé d'un mouvement rectiligne uniforme de 
— 
vitesse uniforme V. 
On a immediatement C, N C,-=@ etune droite d'univers a tous ses points contemporains 
2 2 


ou aucun. 
Considérons maintenant les différentes opérations suivantes : 








7 Peer) IS AE et re = 
y-B 7 x y-x x’ c t y-B-x 
> €. 


Cela montre qu'une droite de C, pour un observateur est une droite de type C, pour tout autre 


observateur non fixe( B > 0) par rapport au premier 


Les 2 sont de type "espace" et pour le deuxième observateur 


; #3 c 
elle represente un mouvement supraluminique V'= — > c. 





y v-B cet 0 0 
(b) | _e& |= re(8- 7) =| cer. Fe: 
rB 7 B B rB 


. F . . C A 
Cela montre qu'une droite de C, représentant un mobile de vitesse V =-— peut etre vue comme 


une droite de type C;. 


Conclusion : 

Identifier les droites affines de l'espace vectoriel réel E de dimension 4 de la relativité restreinte 
aux trajectoires de mouvement rectiligne uniforme est (un peu ) abusif mais vrai si on considère 
plusieurs observateurs non fixes les uns par rapport aux autres et des transformations linéaires : 
Une droite d'univers est vue comme un mouvement rectiligne uniforme par au moins l'un des 
observateurs. 
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